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Abstract: 

This paper shows that the approach based on Gribov's ideas for the photon-proton in- 
teraction, is able to describe the experimental data over a wide range of photon virtualities 

= — 100 GeV^^ and energies ^/s = W = ?> — 300 GeV^. A simple model is suggested which 
provides a quantitative way of describing the matching between short and long distances ( 
between "soft" and "hard" processes ) in photon-proton interaction at high energy. The main 
results of our analysis are: (i) the values of the separation parameters differentiating the "soft" 
and "hard" interactions are determined; (ii) the additive quark model can be used to calculate 
the "soft" contribution to the photon - proton interaction; (iii) a good description of the ratio 
c"l/c"t is obtained ; and (iv) a considerable contribution of the "soft" process at large Q^, and 
of "hard" processes at small was found. 
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1 Introduction 



The rich and high precision data on deep inelastic ep scattering at HERA 0, covering 
both low and high regions, lead to a theoretical problem of matching the non-perturbative 
("soft") and perturbative ("hard ") QCD domains. This challenging problem has been under 
close investigation over the past two decades, starting from the pioneering paper of Gribov ||^ 
( see also [|] ). 

Based on Gribov's general approach, one can interpret two time sequence of the 7*p inter- 
action (see Fig. |I]): 



1. First the 7* fluctuates into a hadron system (quark- ant iquark pair to the lowest order) 
well before the interaction with the target; 

2. Then the converted quark- ant iquark pair (or hadron system) interacts with the target. 
These two stages are expressed explicitly in the double dispersion relation suggested in Ref. 
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where M and M' are the invariant masses of the incoming and outgoing quark- ant iquark pairs, 
cr(M^, M'^, s) is the cross section of a qq interaction with the target, and the vertices r^(M^) 
and r2(M'2) are given by T'^{M'^) = R{M'^), where /^(M^) is the ratio: 

cr(e+e" 1^ 

which has beem measured experimentally. For large masses (M, M' >> Mq) we have r(M^ > 
M^) X r(M'2 > M^) — y R^W^ > M^) = 2, where we assume that that the number 
of colours Nc = 3. Mq is a typical mass (a separation parameter), which is of the order of 
Mo = IGeV, determined directly from the experimental data 

The key problem in all approaches utilizing Eq. is the description of the cross section 
(t(M^, M'^, s). In this paper we follow the approach suggested in Ref. [^], which is based on 
the ideas of Badelek and Kwiecinski [^]. Below we list the mains steps of this approach. 



1. We introduce Mq ^ IGeV in the integrals over M and M' which plays the role of a 
separation parameter. For M, M' > Mq the quark - antiquark pair are produced at short 
distances ( r_L oc < ), while for M, M' < Mq the distance between quark and 
antiquark is too long (rj_ oc -g > jj^^), and we cannot treat this qq - pair in pQCD. 
Actually, we cannot even describe the produced hadron state as a gg - pair; 



1 




a(M,M',s) 



Figure 1: The generalized Gribov's formula. 



2. For M, M' < Mq we use the Additive Quark Model g] in which 

The above assumption allows us to simplify the Gribov formula of Eq. (|lD: 

R{M'^)M'^dM'^ 



37r J (g2 + 



-fTjv(M^s) ; 



(3) 



(4) 



3. For M, M' > Mq we consider the system with mass M and/or M' as a short distance 
quark - antiquark pair, and describe its interaction with the target in pQCD (see Fig. H). 

The exact formulae for a^'""^{M'^, M'^, s) is given below (see also Ref . ) , but one can see 
from Fig. |l|, that this interaction can be expressed through the gluon structure function, 
and it is not diagonal with respect to the masses, contrary to the "soft" interaction of a 
hadron system with small mass. 

4. In principle, the short distance between the quark and antiquark leads to short distances 
in the gluon-nucleon interaction. However, the typical distance of the gluon interaction 
tg oc (see Figs 1 and 2) is larger than the size of the quark - antiquark pair and 
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l± < k±. It turns out that the calculation with Mq ~ 1 GeV, demands a new scale 
in the gluon-nucleon interaction. We introduce it, assuming that the gluon structure 

function xG{x,l'^) behaves as xG{x,ij?) x ^ for Z^. < /^^- This means that we assume 
that the gluon-hadron total cross section is not equal to zero for long distances (in the 
"soft" kinematic region). It should be stressed that we introduce two scales for soft 
nonperturbative interactions, distinguishing between quark and gluon interactions. The 
two scales appear naturally in our models for "soft", nonperturbative interaction. For 
example, in the AQM which we use for 0"*^°^*, these two scales are the size of a hadron 
(distance between the quark and antiquark), and the size of the constituent quark which 
is related to gluon interaction scale. 

5. As has been discussed, we use the AQM [see Eq. (|)] to calculate cr^°^*. For the vector 
meson cross section the AQM leads to 

c^R =^((^(7r+p) + (T(7r>)) . (5) 

For the pion-proton cross section, we use the Donnachie - Landshoff Reggeon parameter- 
ization 1^, with an energy variable which is appropriate for the interaction of a hadronic 
system of mass M with the target (see for details). 

aM^N{s) = A( — )"-W-i + B (^)"«(o)-i, (6) 
Xm Xm 

with 



ttjp = 


1.079 


OllR = 


0.55 






Xm = 






s 


A = 


13.1 mb 


B = 


41.08 mb 



(7) 

The values of A and B are chosen so that (|]) is valid for the p-proton interaction. Using 
d^) in Gribov formula (^) we derive the soft transverse contribution to cr(7*j9): 



3tt Jo (Q2 + M^y I Xm xm 



It has been shown in that even though the simple model discussed above reproduces the 
main features of the experimental data, it has several deficiencies which we wish to address: 
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Figure 2: "Hard" contribution to a{M'^ , M''^ , s) in perturbative QCD . 



• The calculations of the "soft" cross section appeared to overestimate the data, and to 
overcome this, the soft contribution was multiplied by a factor k, = 0.6. There is no 
physical justification for such a small value. 

• The contribution of a longitudinal polarized virtual photon to cr{'y*p) was neglected. 

• An old GRV'94 parameterization was used for the gluons distribution inside the proton. 
This gave an energy dependence of <j{Yp) which is steeper than the published data. 



In the present paper we reexamine these points developing a formalism which also takes 
into account ul, the longitudinal part of a{'-f*p), for both the soft and the hard components. 
We find that the contribution of ul significantly improves the energy dependence of (t{'J*p), in 
a good agreement with the data for x < 10~^. 



A similar approach to photon - proton interaction was also developed in Ref.|TO[ where 
quite different scales have been introduced. The goal of such an approach (see Refs.[^],0 
and |TD|) is to find parameters that separate the long distance (non-perturbative) and short 



distance (perturbative) interactions in QCD. A similar philosophy to ours is used in Ref.|jTT 
where an attempt was made to describe photon - proton interactions, assuming that the main 
contribution stems from short distances, and the long distance interaction appears as a result 
of the shadowing corrections. 

The description of criY'p) can be achieved in quite a different way using the generalized 
vector dominance model @] or a Regge motivated fit of the experimental data[|l2] ||l3l. In these 



approaches one looses the explicit connection with the microscopic theory, which makes futher 
theoretical interpretation of the results rather difficult. 
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Figure 3: Diagrams which cntribute to the scattering of a qq pair off the target proton. 



2 Description of the Model 



2.1 The Contribution of a Transverse Polarized Photon 

As we have discussed previously, the main assumption of our mo del is that a(M2,M'2,s) m 
Eq. (|1|) can be calculated as 

a(M2, s) = a'j^f\M^ < s) 6{M^ - M'^) + a^'^'^M^ > M'^ > s) . (9) 

In the previous section we have discussed the calculation of 0"''°^*. Now we present, for com- 
pleteness, our formulae for a^^^'^ which have been given in Ref . . The pQCD contribution for 
o"T, within the framework of the two gluon model, is illustrated in the four diagrams shown in 
Fig. |[ The production amplitude Aiw of the qq, can be factorized into the wave function ipxy 
of the qq system inside the virtual photon, and the amplitude Txy for the scattering of the qq 
pair off the proton: 



Jo 



Mxx'ik^,z) = ^N, 

The transition ampiltude Txy for the four contributions is: 
Txyik^', k±) = 

r rfifi . r 

2(5(fcl - ki_) - 5{k'^ -k±- i±) - 5{k\ 



(10) 
(11) 



Arrs r (P 



2Nr 



as(£i)/(x,£l) 
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where, 



f{x,f 



dxG{x, 
(9 In £2 



see 



H| for details). Substituting (^Tj) in (|T0|) we have 



(12) 



(13) 



where the definition of Aip is to be understood from Eqs. (|10|)-([T3|): 

= 2ip{k±, z) - tp{k^ - i^, z) - i){kx. + z) . (14) 

'j/'AA' has been calculated in both for a transverse and for a longitudinal polarised incoming 
photon. The wave function for a transverse polarized photon and the amplitude for producing 
qq with helicities A and A' read: 

2e± ■ 



V^AA' (k±,2) 



-A^AA' 



A,-A 



Z/e [(1 - 2z)\ T 1] 



2 ' 



(15) 



-2Zfe U-j=] 5a,-a' [(1 - 2z)\ T 1] e± ■ ^«s(^l)/(x, ^1) x 
1 1 Q^-k^^ + il 



2k, 



2k I 



g2 + kj_^ + 



(16) 



In ([TSD and (|T6D, Zj is the charge of the quark with fiavour / in units of the electron charge — e, 
Q"^ = z{l — z)Q'^ and e± denotes the photon polarization vector presented in a circular basis, 



e± = ^(0,0,l,±i). 



(17) 



To evaluate the cross section one should, first, sum over the quark helicities A and A', and 
average over the two transverse polarisation states (±) of the photon 



\M 



T\2 



327rV 



z^ + il-zf 



k^' - 



+ 



Q^-k^' + il 



> . 



The cross section is obtained by integration over z and over 



A'k'^Nc Jo 



dz 



z' + {l- zf 



k^' - 



Q^ + k] 



de 



+ 



Q^-k^' + il 



+ kj_^ + 



ik 



2,^2 



(19) 



Following , in order to perform the integration over we introduce the variables M and M: 



M 



2 



z{l - z) 

= — (20) 
zil - z) 



and rewrite the integrals of ctt in terms the variables M, M and 

dM^ rdM^ rdB 1-2 



2 



■2 

1-4^ 



M^-g^ ^ g2_^2^^2 ^^^^ 

+ '^^ y(Q2 + M2 + M2)^ - 4M2M2 

Using a generalization of the Gribov approach, we take a cutoff on the integration over M2 and 
insert the ratio @ inside the integration sign, 

.hard _ «em dM' R{M') f dM' f dil 1-2^ ^ ^^^^^^^ 



(22) 



M2 

M2 - Q2 ^ Q2_m2^ ^2 

+ y(g2 + M2 + M2)' - 4M2M2 



The last integral in ( p2D can be integrated by parts, using (|1^). In the limit 4i\ ^ M'^, the 
value of the integral is dominated by the upper limit of the integration, therefore we replace 
in as(^i) and xG^Xji^) with M'^/i, and obtain, for three colors 

h,rd 27iaem rdM^ RiW) foo dM^ M\ M\ 



+ M'^ Jo 4 ' ' ' 4 

M2 - g2 g2 ^ jQ-2 _ j\^2 



M2 + g2 J(g2 + M2 + M2)2 - 4M2M2 



• (23) 

In the case of heavy quarks, we assume that the quark is heavy enough so that any contribution 
to the soft cross-section can be neglected. The heavy quarks contribution is then written with 
the replacements: 4M2 M2 _^ 4 (j\^2 _ 4^^^ j^2^ g^^^ R{Q'^) R^^iW^). R^^ is the heavy 
quark contribution to the ratio (0) and mg is the mass of the heavy quark, 

<^T,QQ = -^r- ... r^,. , L ^>Z"s(^)a;G(x,— ) 



^'^^ 3 A/2 g2 + M2 io M4 



+ ^( g2 + M2 + M2 )2 _ 4 (M2 - Ami) 

In the above formulae x = x(M^) = (Q^ + M^)/14^^ with being the center of mass energy 
of the photon-nucleon interaction. 

2.2 The Contribution of a Longitudinal Polarized Photon 
2.2.1 "Soft" Contribution to (x^ 

A priori, it is straight forward to write the formula for the "soft" component of ctl in AQM. The 
result is similar to our model for a^, except that in the numerator should be replaced with 
the photon virtuality. This factor causes the longitudinal cross section to vanish for — > 
(see d^) for the values of the parameters). 

L 3vr Jo (Q2 + JVf2)2 I ^Xm^ ^Xm' ) ^ ' 

However, it turns out that the AQM over estimates the experimental data, and it has 
to be reduced with some phenomenological procedure, such as a numerical factor or a 
phenomenological function which decreases with [IC]. In our formalism the ratio between 
the AQM and pQCD contributions depends on the separation parameter Mq. Thus, by lowering 
the value of Mq the "soft" component is suppressed. We have found that for the longitudinal 
part of the cross section, taking any value of Mq which is below the resonances mass fits the 
experimental data, as opposed to the transverse cross section where we used 0.7 < Mg < 
0.9 GeV^ . 



2.2.2 "Hard" Contribution to 



We calculate the pQCD contribution for ctl, from the diagrams of Fig. 0, using Eqs. (p!OD-(p^ 



For the case of a longitudinal polarized incoming photon, we use the expression derived in ||T4 
for the wave function 

1 



-25 



ZfeQz{l — z) 



(26) 



and substitute it in 
helicities A and A' 



to obtain the amplitude for a longitudinal photon to produce qq with 



M 



L 

AA' 



-46x,-x'ZfeQz{l 



z]i- 



1 



g2 + kj_^ + 



4k 



(27) 
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The cross section is then obtained by squaring, summing over the final hehcities and integrating 
over z and k \ , 



Nr. 



■-asiil)fix,il)x 



(g2 + fcj 



2 _L ^2 r _ 4A;_L^^2 



(28) 



Changing the integration variables fc^. ^1 M and M, respectively, using (^), we have 



AT. 



1 



Q2 + M2 .//g2 + ]Vf2 + M2^2_4]Vf2]Vf2 



(29) 



We now write the formula for o^^'^^^ in the same way that we did in section by taking a 
cutoff on the integration over M^, and insert the ratio (^ inside the integration sign. 



„hard 



Ml Q^ + M^ Jo M4 Jo 
1 



dzf (x, z{l - z)M^) as {z{l - z)M^) x 
1 



Recalling (|12|) , we can integrate Eq. ( ^0|) by parts and obtain (for A/'c = 3), 

R{M^)dM^ f^dM^ 



(30) 



^hard 



1 1 



xG (x, M^) 



Q"^ + Y^(g2 + M2 + M2)2 - 4M2M2 

M2 (^Q2 ^J^2 



(g2 + M2 + M2)2 - 4M2M^ 



(31) 



where. 



xG (x, M^) = f xG (x, z(l - z)M2) as (;z(l - 2)M2) dz . 



(32) 



For <y^QQ, the contribution of heavy quarks to the longitudinal component of cr(7*p), we replace 
iM'^M'^ by 4{M^ - Am^Q)M'^ and R{Q^) with R'^'^{M^), as we did for the transverse part: 



a 



hard 
h,QQ 



cm ^2 



1 



1 



Q^ + M^ ^^Q2 + M2 + M2)2 - 4(M2 - 4m^)M2 



M2 (Q2 + - 



(g2 + M2 + M2)2 - 4(M2 - 4m|)M^ 



(33) 



3 Comparison With The Experimental Data 

We now present the results for our calculation of cr(7*p) from the master equation 

^(^,*^\ ^soft I „hard i „hard i „soft i „hard i „hard {'3A\ 

As stated, our model has three parameters, namely, Mq^l, Mq^t and /x. In [|| the relatively 
high value of Mq = \/5 GeV was taken, in order to give an energy dependence that was 
in a reasonable agreement with the published experimental results. We found that if the 
hard longitudinal contributions are not neglected, the transverse separation parameter can be 
reduced to the value of 0.7 < Mq^ < 0.9 GeV^, while the longitudinal one can be taken to be 
Mq ^ 0.4 GeV^. The results of the calculation are stable if the separation parameters are 
chosen within these bounds. 

In the calculation of the hard components of a{'~f*p), we need to specify the input gluon 
distribution xG{x, Q"^) which appears in the formulae. We have tried several options: GRV'94 
rg, GRV'98 |T3 and MRST'98 We compare the results of each input distribution for 



the calculation of cr(7p) in Fig. ^, where it is obvious that the best parameterization for our 
purpose is MRST'98 which yields, together with our formalism, a good description of the data. 

A difficulty in our calculation comes from the integration at very low M^, where the pub- 
lished parameterizations for xG are not valid. Naturally one has two possible options: the first 
is to impose a low cutoff on the integration below = /x^. The second option is to use the 
general property of the gluon structure function which is linear in at very small values of 
g^, and to rely on this property for the low mass region: 

xG(x,l'^ < fi'^) = —xG(x,fi'^) . (35) 

/i2 
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Figure 4: Comparison of the gluon distribution parametrizaions which were used in Eq. (34) 
to calculate a{'y*p). The full triangles correspond to data points which have been 
extracted in ^l^J from experimental data. 
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In our calculation for the low mass integration we used (|35| ) for the gluon distribution. The 
coupling as(Q^) was also kept fixed for < /x^. The value of yU^ was taken as the minimal value 
which is allowed in the input parameterization that had been used: fi^ = 0.4, 0.8, 1.25 GeV^ 
for GRV'94, GRV'98 and MRST'98 respectively. Notice that for the MRST parameterization 
jj? = 1.25 GeV^ > Mq. This means that the transition from soft to hard is not sharp, and we 
still have ( pSf ) as a soft signature inside our hard formalism at Mq < < /i^. 

In Fig. 1^ we show the results of our calculation for a{'~f*p) as a function of W'^ for fixed values 
of Q^, together with the published experimental results. The calculated results are shown as a 
band which corresponds to the two limits on Mq.t- It is clear from the figure that the width of 
the band decreases with increasing Q^. As for the longitudinal separation parameter, the results 
are not sensitive to the choice of Mq^l, and in all our calculations we simply used Mq^l = 2m^. 
The thick line marks the region of small x: to the left of it, x is not small enough to justify 
using our model. One can see that for x < 10"^, our results reproduce the experimental results 
both in value and in energy dependence. 

In Fig. P we show our calculations for o"(7*p) for fixed values of W"^, as a function of Q^, 
scaled by factors of 1 — 128. The small vertical thick lines delineate the boundary x = 10~^, to 
the right of the marks our results are not reliable. 

The longitudinal component of cr{'y*p), according to our calculation, reaches a maximum 
value of 25% of the total cross section at x ~ 10~^. For small x, this result is closer to the 
published data 120] than the results of Ref. [10], where they obtained ctl for x ~ 10~^ to be 



15% of the total cross section. The ratio of the longitudinal cross section to the transverse cross 
section, as a function of x, is shown in Fig. ^ together with the (small x) data of Ref. |2^ . 



Fig. H illustrates the importance of the heavy quarks contributions o'^^^qq and ct^^^q [see Eqs. 
4]) and (0)]. The contribution has a mild W^-dependence but it is Q^-dependent, and for 



large Q^, the ratio ci^'^^^^/o"''^'^* gets as high as 0.3 — 0.4 . Another ratio which we present is the 
ratio, of the hard to the soft contributions. Since Mo,l is smaller than the lightest resonance 
mass, the AQM contribution to the longitudinal component is supressed, and therefore what 
we call "soft" contribution is actually the transverse AQM calculation. 

In Fig. ^ we present the ratio o-soft/o-hard fixed values of Q^, where a clear power law 
behaviour as a function of W"^ can be seen. In soft processes, the hard contribution is minor, 
but is still present even for relatively small values of Q^. The soft signature in hard processes is 
also present, but is smaller. For high virtualities it decreases from a few precent at low energy, 
to less than 1% at high energies. 



4 Conclusions 

In this paper, we have shown that an approach based on Gribov's proposal, provides a suc- 
cessful description of the experimental data on photon-proton interaction, over a wide range 
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15 



0.5 



Q2 = 65 
q2 = 35 



Q2=15 



q^ = 8.5 
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of photon virtualities < < 100 GeV^ and energies 3 < ^ {= W) < 300 GeV. The key 
assumption on which our aproach is based, is that the non-perturbative and the perturbative 
QCD conributions in the Gribov formula can be separated by the parameter Mq. Oue successful 
reproduction of the experimental data (see Figs. ^ and shows that this assumption appears 
to be valid. It further lends credence to using the additive quark model (AQM) to describe 
the non-perturbative contribution. The successful use of the AQM leads to an improved result 
presented in this paper, compared to our previous result 0], where we found it necessary to 
introduce an damping factor for the AQM. 

The second important by-product of our calculation is the simple method we have used 
to separate between non-perturbative ("soft") and perturbative ("hard") contributions. We 
have used three separation parameters: Mo,t, ^o,l and /i. The values of these parameters, 
= 0.7 - 0.9 GeV^ M^^ < 0.4 GeV^ and /x^' ~ i GeV^ were determined by fitting to 
the experimental data. We believe that these values may be useful in the future, for more 
theoretical description of the matching between "soft" and "hard" processes in QCD. 

The third result which we find interesting, is that the GRV parameterizations of the structure 
functions cannnot successfully describe the experimental energy dependence of the total cross 
section at small values of (see Fig. while the MRST parameterizations can. This result 
shows the interdependence of deep inelastic scattering data, and the theoretical description of 
the mathching of the "soft" and "hard" contributions. 

In addition we obtain: 

1. A description of the ratio ctl/ctt (see Fig. |^) which is in good agreement with the ex- 
perimental data and other approaches. The fact that appears to have only a "hard" 
contribution should be stressed, as this could be a possible window for particular features 
of non-pertutbative QCD; 

2. A considerable contribution of the "soft" processes at rather large values of photon vir- 
tualities For example at = 15 GeV^ a'°^/a^''"^ ^ 4% at W = 10 GeV (see 
Fig. I). We also find, at low Q'^, a contamination of the "soft" processes by the "hard" 
ones. For high energy real photoproduction this contamination amounts to about 5% . 
We believe that this fact should be taken into account when interpreting experimental 
data, especially as far as their energy dependance is concerned; 

3. The prediction of the cross section for the heavy quark production (see Fig. ^ which is 
almost W independent, while showing a steep decrease. 

We propose a simple and successful phenomenological model for the photon-proton interac- 
tion, which provides a method for matching the "soft" and "hard" interactions at high energy, 
and which can be a guide for future theoretical approaches. It is worthwhile mentioning, that 
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even in the present form our model can be useful for determining the initial parton distributions 
at leading twist, for the DGLAP evolution equations. 
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